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ABSTRACT 

We studied the effects of rotation on standing accretion shock instability 
(SASI) by performing three-dimensional hydrodynamics simulations. Taking into 
account a realistic equation of state and neutrino heating/ cooling, we prepared a 
spherically symmmetric and steady accretion flow through a standing shock wave 
onto a proto-neutron star (PNS). When the SASI entered the nonlinear phase, we 
imposed uniform rotation on the flow advecting from the outer boundary of the 
iron core, whose specific angular momentum was assumed to agree with recent 
stellar evolution models. Using spherical harmonics in space and Fourier de- 
compositions in time, we performed mode analysis of the nonspherical deformed 
shock wave to observe rotational effects on the SASI in the nonlinear phase. We 
found that rotation imposed on the axisymmetric SASI did not make any spiral 
modes and hardly affected sloshing modes, except for steady / = 2, m = modes. 
In contrast, rotation imposed on the non- axisymmetric flow increased the ampli- 
tude of spiral modes so that some spiral flows accreting on the PNS were more 
clearly formed inside the shock wave than without rotation. The amplitudes of 
spiral modes increased significantly with rotation in the progressive direction. 

Subject headings: hydrodynamics — instability — neutrinos — supernovae: 
general 
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Introduction 



Core-collapse supernovae are among the most energetic explosions in the universe, cat as- 
trophically destroying massive stars. Because they are relevant to many astrophysical phe- 
nomena (e.g., formations of compact stars like neutron stars or black holes, nucleosynthe- 
sis, neutrino and gravitational emissions), their physics has been of wide interest to the 
astrophysical community. Regardless of rigorous studies on core-collapse supernovae, the 
explosion mechanism is still not completely understood. Except for the lower mass progen- 
itors, spherically symmetric supernova simulat ions have not yet produced explosions (e.g., 
Liebendorfer et al.ll2005l ; ISumiyoshi et al.ll2005l and references therein). As a result of cur 



rent observations revealing the aspher ical nature of the explosion (e.g., IWang et al.l 12002 



Tanaka et al.ll2007l ; iMaeda et al.ll2008l ). multi-dimensional studies and simulations of core- 
collapse superno vae have explored various mechanisms o f asphericity, such as t he roles of 



convection (e.g.. iHerant et al.lll994t IBurrows et al.lll995l ; iJanka fc Mueller! Il996l ). magnetic 



field and rapid rotation ( e.g.. iKotake et al. 



shock instability (SASI) ( Blondin et al 



20061 : lOhnishi et al.l [2009 . 
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20081 ). and g-mode oscilla- 



SASI, the hydrodynamic instabilit y of a standing shock wave, was origin ally studied 
in the context of accreting black holes (iFoglizzol l200ll ; iHouck fc Chevalier! If 9921 ) . The main 
characteristic of SASI is that lower I modes dominate the flow dynamics, where I stands for 
the polar index of the spherical h armonics Yj m . The im portance of SASI to the supernova 
dynamics was first pointed out by lBlondin et al. (b()03h . who demonstrated that non-radial 
perturbations added to a standing supernova shock wave grow exponentially with time and 
lead to I — 1 or 2 mode deformations in the linear phase, and then the slos hing motion of 



the sh ock wave induced more violent turbulent flows in the nonlinear phase. lOhnishi et al. 



( 120061 ) indicated that SASI may be the key to the explosion mechanism because SASI can 
decrease critical neutrino luminosity for the shock revival. To date, most realistic axisym- 
metric two-dimensional core-colla pse simulations have revealed the appearan ce of SASI in 
the initial phase of the explosion (IBurrows et al.ll2006l ; iMarek &: Jankall2007l ). Moreover, it 
has been suggested that such a lower mode explosio n is favorable for reproducing the syn- 
thesized elements o f SN19 87A (Kifonidis et al.l 120061) an d for explaining the origins of kick 



(jScheck et al.ll2004l) . spins (IBlondin &: Mezzacappa 
2008|) of pulsars. 



20071 ). and magnetic fields (jEndeve et al. 



Reflecting the growing importance of SASI, the physics behind SASI in supernova cores 
has recently drawn much attention. Two types of mechanisms are considered: the advective- 
acoustic cycle (AAC) and the purely acoustic cycle (PAC). In the AAC scenario, the entropy 
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and vorticity perturbations advected toward PNS generate a sound wave at the location of 
the largest velocity gradient of the stationary flow. The sound wave propagates toward the 
shock wave and distorts its configuration, depending on the non-radial distribution of the 
fluctuation pressure. The de formed shock wave i nduces further amplification of the entropy 
and vorticity perturbations (jFoglizzo et al.l 120071 ). For the PAC scenario, the standing pres- 
sure wave propagates in the circumferential direction in the region between the spherical 
accretion shock wave and PNS. When the postshock pressure is slightly higher than unper- 
turbed pressure, it pushes the shock wave outward. The outward displacement of the shock 
wave leads to an increase of postshock pressure in the inner region, while the postshock 
pressure immediately behind the shock wave decreases. Th us, the amplitude of the pressure 
fluctuation increases further (IBlondin &: Mezzacappall2006l ). At this time, which mechanism 
really works in the supernova cores is still a topic of debate. 

For the past few years, the no n-axisymmetric features of SASI have been invest igated 
with three-dimensional simulations (IBlondin fc Mezzacappdl2007l ; llwakami et al.ll2008l ). The 
modes of SASI are divided into sloshing modes and spiral modes. Axisymmetric m = 
modes thus far studied in two-dimensional axisymmetric models (symmetric axis is z-axis) 
and degenerate \m\ = I modes (symmetric axes are x-axis, y-axis, and so on) are classified 
as sloshing modes, where m stands for the azimuthal index of the spherical harmonics Y™. 
When we impose random perturbation or rotating flow on these axisymmetric flows, the 
degeneracy is broken and the rotational modes emerge. In this situation, the +m modes have 
different amplitudes than the — m modes. Such rotating non-axisymmetric m^fl modes are 
called spiral modes. The sloshing modes can be expressed as degenerate spiral ±m modes 
which have the same amplitudes as each other. The growth of spiral modes in the linear phase 
has been examined by sim ulation with a two-dim ensional polar grid of a thin wedge over 
the entire equatorial plane (IBlondin fc Shawl 120071 ). A deformed shock wave resulting from 
the growth of the spiral modes produces a spiral accretion flow inside the shock wave, which 
transfers the angular momen tum onto the PNS. A three-dimensional SASI simulation by 
Blondin fc Mezzacappal (120071 ) demonstrated that the spiral mode of m = 1 grew dominantly 
and generated a strong rotational flow, continuing to develop in the nonlinear phase when 
random perturbati ons were imposed on the non-rotating or rotating progenitor. However, 
our previous work ( Iwakami et al.l2008 ) revealed that the rotational flow did not develop as 
much for non-rotational models. An equi-partition was nearly established among different 
m modes on a time average in the nonlinear phase. In the flow, high-entropy blobs were 
produced inside the shock wave. High-velocity accreting matter ran outside the blobs, and 
circulating flows formed inside the blobs. The high-entropy blobs repeatedly came and wen t 



during the nonlinear phase. In contrast with the outcomes of IBlondin fc Mezzacappal (2007). 



our results were similar to those obtained by the Garching group (jWoosley &: Janka 



2005|), 
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which revealed that large bubbles of radiation formed inside the shock wave. Why these 
results differed so much is not yet understood. 



The SASI for the rotating progenitor model has also been studied. iBlondin &: Mezzacappa 



(120071 ) demonstrated that the rotation of the infalling gas helped to initiate the spiral modes 
of SASI, and t he rotati o n axi s of the flow was roughly aligne d with the spin axis of the 
progenitor sta r. lLamind ( 20071 ') and lYamasaki fc Foglizzol (120081 ) analyzed the effect of rota- 
tion on SASI. lLamineJ toO^ i deriv ed an approximate dispersio n relation for oscillations of 
the spherical accretion shock wave. lYamasaki &: Foglizzol (120081 ) used perturbative analysis 
with WKB approximation for the cylindrical shock wave. Although they devided over the 
opinion which AAC or PAC worked dominantly, they found that the growth rate of the 
modes rotating in the same direction as the flow was increased by rotation. In the present 
study, we investigated the effect of rotation on SASI in the nonlinear phase. We introduced 
a rotation whose axis was along the z-axis into the perturbed flow in the nonlinear phase. 
We conducted mode analysis expanding the deformation of the shock wave with spherical 
harmonics in space and Fourier series in time to distinguish between +m and — m modes. If 
the fe ature s in the nonlinear phase we re similar to the results of linear analysis by [Laming 
(120071 ) and lYamasaki fc Foglizzol (120081 ). it was expected that for clockwise rotation of uj^ < 
the m < modes should grow much more than the m > modes. 

The organization of this paper is as follows, we describe the numerical models and the 
formulations for mode analysis in §2, we present the results of computations in §3, and we 
make a summary and discussion of this study in §4. 



2. Methods of Computation and Analysis 



2.1. Numerical Models 



The numerical method was exactly t he same as that use d in our previous paper (llwakami et al. 



20081 ). Employing the ZEUP-MP code (IHayes et al.ll2006l ) for the hydro solver, we solved 
the dynamics of an accretion flow of matter attracted by the PNS and irradiated by neu- 
trinos emitted from t he PNS. The Shen EO S (jShen et al.lll998l ) was implemented according 
to the prescription in iKotake et al.l (120031) . To treat the ne utrino heating and cooling, we 
used the light bulb approximation (see lOhnishi et al.l 12006k for details), in which the neu- 
trino heating is estimated under the assumption that neutrinos are emitted isotropically 
from the central object with a fixed neutrino flux, and that the matter outside PNS is opti- 
cally thin. For simplicity, we considered only the interactions of electron-type neutrinos and 
anti-neutrinos. Their temperatures were constantly assumed to be the typical values in the 
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postbounce phase (i.e., T Vc = 4 MeV and T Pc = 5 MeV). The neutrino luminosity was fixed 
at L v = 6.0 x 10 52 ergs s _1 . The mass of the central objects was assumed to be M in = 1.4M , 
and the mass accretion rate was set to be M = 1 M s _1 . 

A staggered mesh in the spherical polar coordinates system was used. The mesh had 
300 radial mesh points to cover r- m < r < r out , and 30 polar and 60 azimuthal mesh points 
to cover the entire solid angle. Here, r out = 2000km was the radius of the outer boundary 
at which the flow was supersonic, and r in ~ 50km was the radius of the inner boundary 
located roughly at the neutrino sphere. This angular resolution w as enough to investi gate 



the characteristics of SASI being dominated by lower modes (see llwakami et al.l 120081 . for 
resolution tests). 

For the boundary conditions, we imposed the fixed-inflow condition on the outer bound- 
ary, and the free-outflow condition except for the radial velocity v r on the inner boundary. 
At the points on the inner boundary, v r was fixed with the value of initial flow in order to 
obtain a steady flow at the onset of calculation. At the points on the inner ghost boundary, 
v r was determined to satisfy v r>Q = v T ^r\jr\ (v r>0 > iv,i), which meant the conservation of 
mass flux for ID flow, where v r $ is the v r at r and v r> i is the v r at r\ (r < ri). We con- 
firmed the influence of inner ghost boundary condition to change the relational expression to 
the fixed-outflow condition (tv,o = u r ,i) or the fixed-initial-flow condition (tv j0 < v r ,\)- The 
change of relational expression in the ghost boundary made little difference in the essential 
features of the flow. 



We used the spherically symmetric steady flow as the initial condition (jYamasaki fc Yamada 



20051 ). The radial di stributions of various variables for the unperturbed flows were given in 



our previous paper (llwakami et al.l 120081 ) . In order to induce non-spherical instability, we 
added radial velocity perturbation, 5v r (9,(p), to the steady spherically symmetric flow ac- 
cording to the following equation: 

v r (r,9,<!>)=v 1 r D (r)(l + 6v r (9,<t>)), (1) 

where vl D {r) is the unperturbed radial velocity. In this study, we considered two types of 
perturbation: (1) an axisymmetric I — l,m — single-mode perturbation, 

ry 

5v r (9, 0) cx \ — cos9, (2) 

V 47T 

and (2) a non-axisymmetric random multimode perturbation, 

Sv r (9, (f>) oc rand (0 < rand < 1), (3) 



where rand is a pseudo-random number. These perturbation amplitudes were set to be less 
than 1% of the unperturbed velocity. 
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2.2. Introduction of Rotation 

In contrast with a spherically symmetric accretion flow, the method to construct the 
stationary rotational accretion flow running through the rotationally deformed shock for 
core-collapse supernovae has apparently not been developed thus far. Hence, to investigate 
the effects of rotation, we imposed rigid rotation on only the outer boundary. The imposed 
rotation is advected toward the shock wave by the accreting flow. We focused on only the 
nonlinear phase in this study; therefore, we introduced a rotational flow into the flow field in 
which SASI had already developed in the nonlinear phase. We then analyzed the difference 
between the result with the introduction of the rotation and that without. 

We yielded the rotation described as follows on the outer boundary: 

vl D (r,8)=vl D (r)(3^m8, (4) 

where v^ D (r, 6) denotes the unperturbed component of velocity and denotes the rotation 
parameter. We examined the flow characteristics for (3$ = 0.005 — 0.050, which corresponded 
to the specific angular momentum L ~ 0.2 — 2.3 x 10 16 cm 2 s -1 on the equatorial plane. 
Table [T] lists all the models used in this study. The L of these models was chosen to be 
reconciled with th e results of t he re cent presupernova calculations of rotating stars with 



magnetic fields by lHeger et al.l (120051 ). 



Figure [T] illustrates the evolution of radial distributions of along the equatorial plane 
{9 = 90°, = 0°) for Model Al, by which one can see how the rotation given to the outer 
boundary was advected toward the PNS. At first, was zero everywhere when the rotational 
flow for Pfj, = 0.01 was introduced into the outer boundary at t — 0ms. The rotational flow 
with v<t, was then advected from the outer boundary to PNS with increasing its value by 
compression for t = 20 — 100ms. At t — 100ms, the rotational flow reached the surface 
of PNS and hardly changed for t = 150 — 300ms, where the plots for t = 200ms and 
t = 300ms overlapped. 

Figure [2] depicts the time evolutions of the average shock radius Rs for unperturbed 
flow. The average shock radii for Models AO, Al, A2, and A3 were almost constant from 
150ms to 300ms. The average shock radii for the rotating flow of Models Al, A2, and A3 
were slightly larger than that for non-rotational Model AO because of the centrifugal force. 
When we applied more rapid rotation like Model A4, the shock wave continued to expand 
outward. The models for (3$ < 0.015 were the non-exploding models for unperturbed flow, 
and allowed us to analyze clearly the role of rotational SASI in flow dynamics. Thus, we 
focused mainly on the models for /3$ < 0.015. 

Figure [3] displays the side views of the iso-entropy surfaces with the velocity vectors in 
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the equatorial section. The entropy is expressed in the eight iso-surfaces shifting from white 
to the red end of the spectrum with increasing value, and the velocity is denoted by the 
elongated trigonal pyramids shifting from white to blue. The hemispheres (tt/2 < 9 < n) 
of eight iso-entropy surfaces are superimposed on one another, and the outermost surface 
almost corresponds to the shock front. The central region represents the physical quantities 
on the spherical surface corresponding to the inner boundary. An unperturbed flow remained 
spherically symmetric at least until t ~ 300ms (Fig. [3] (a)). When we introduced the clockwise 
rotation into the unperturbed flow at t — 0ms, the rotating flow was advected inward and 
arrived at the inner boundary at t ~ 100ms. Following that, even though the shock wave 
slightly oscillated, the flow remained axisymmetric at least until t ~ 300ms. This result 
indicated that no spiral modes were generated numerically (Fig. [3] (b)). 



2.3. Formulation for Mode Analysis 

This section explains the mode analysis method used in this study. The deformation 
of the shock surface could be expanded as a linear combination of the spherical harmonics 
components F z m (^,0): 

oo I 

Rs(e,M = J2Yl crmriPA), (5) 

Z=0 m=-l 

where Y l m (9, 0) is expressed by the associated Legendre polynomial P; m (cos6 l ) and a constant 
A™ given as 

Y l m (9, 0) = A^P/^cos 9) e im ^, (6) 



2l + l (l-m)\ 

Kl = V 4vr (l + m)V (7) 
The expansion coefficients were obtained by 

/•2-7T pTT 

c?(t)= d<j) d9 sin 9 R s (9, <p,t)Y l m * (9, 0), (8) 
Jo Jo 

where the superscript * denotes complex conjugate. Instead of the expansion coefficients 
c™(t), we used normalized amplitudes c™ (t) / 'cg(t) for analysis. 

Moreover, normalized amplitudes c^(t)/co(i) could be expanded to the Fourier series as 
follows: 

/■oo 

cT(t)/c° (t)= / ^c[»e-^ (9) 
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where u is a real number denoting an oscillation frequency. Thus, we could rewrite (jSJ) as 

oo I 



^ * poo 

R s (9,M/c° (t) = J2J2 / duc?{u)K™P™{cos9)e 

1=0 m=-l J -°° 

Therefore the Fourier expansion coefficients c™(uj) were calculated as 



(10) 



— fdt I d(j) I d9 sin9 [R s (e,(f>,t)/c° {t)] K\ n Pf 1 (cos 9) e *( wt -^), (11) 

— ts Jt s JO JO 



C, (UJ) — 

1 K J U-t 



where t s is the starting time of sampling, and t e is the ending time. The Fourier expansion 
of c™ {t) / c^it) allowed c™(u;) to distinguish between +m and — m modes. 



3. Results 

In this study, we introduced the rotation described above into axisymmetric (2D) and 
non-axisymmetric (3D) flows in the nonlinear phase at t = 400ms to observe the behavior of 
the flow in only the nonlinear phase. We made both 2D and 3D flows in the nonlinear phase 
with three-dimensional simulation as follows. Figure H] presents the plots of the normalized 
amplitudes \cf{t)/cl(t)\ as a function of time for 2D Model Bl and 3D Model CI. In 2D 
Model Bl (Fig. H] (a)), the axisymmetric I — 1, m = perturbation w as imposed on the 



initial flow. As we described in our previous paper (jlwakami et al.ll2008l ). the amplitude of 
the / = l,m = mode grow exponentially in the linear phase, and the growth of amplitudes 
of all modes is saturated in the nonlinear phase. In 3D Model CI (Fig. 0] (b)), the non- 
axisymmetric random perturbation was imposed on the initial flow. As we mentioned in 
the same paper, the amplitudes of all the modes grow from the beginning, and then the 
flow enter the nonlinear phase. For both 2D and 3D flows, we imparted the rotation on the 
outer boundary at t = 400ms (dotted line). At that time, the flow was completely in the 
nonlinear phase. The rotational flow was carried by an accretion flow and arrived at the 
inner boundary around t = 500ms (dashed line). Following that, the flow inside the shock 
wave changed, reflecting the infalling rotational flow. We defined the period from 600ms 
(dot-dashed line) to 1000ms as the rotation phase. We integrated the shock radii Rs(9,(j)) 
from 600ms to 1000ms in Eq. ffTTl) for the Fourier expansion. 



3.1. Axisymmetric Rotational Model 



First, we will discuss the results of axisymmetric flows. All the simulations were carried 
out in three-dimensions even for the axisymmetric flow. We dare to term the axisymmetric 
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flow as the 2D flow not to confuse axisymmetric flow with axisymmetric mode. When we 
introduced the I — l,m — perturbation into the non-rotational flow, only m = modes 
grew, and the flow retained symmetry with respect to the z-axis even i n the nonlinear phase . 



In other words, the sloshing modes grew, but the spiral modes did not (llwakami et al.ll2008l ). 
Confirming that only sloshing modes appeared even for the rotational models, hereafter, we 
focused on the effect of rotation on the sloshing modes in the nonlinear phase. Furthermore, 
we examined that the results of the new mode analysis presented in Section 2.3 to expand 
the geometry of the deformed shock wave with spherical harmonics in space and Fourier 
transform in time, which were consistent with those of the previous mode analysis to expand 
with only spherical harmonics in space at the given time. 

Figure depicts the time evolutions of the average shock radius R$ for 2D Models BO, 
Bl, and B2. The average shock radius initially increased with the growth of SASI, and 
then it remained roughly constant until 500ms with little fluctuation. After 500ms, rotation 
affected the flow dynamics for Models Bl and B2. The magnitudes of average shock radii 
were not as different between non-rotational Model BO and rotational Models Bl and B2. 
We observed that the shock wave burst suddenly if we imposed more rapid rotation on the 
flow. In this paper, we focused on only the non-exploding Models Bl and B2. 

The flow fields for rotational 2D Model Bl are presented in Fig. [6j The velocity vectors 
are superimposed on the cutting plane of the partial cutaway view of the iso-entropy surfaces. 
As with the non-rotation, the flow field retained symmetry with respect to z-axis even in the 
flow rotating globally (Fig. El (a))- No spiral flows formed, so the flow just rotated clockwise 
(Fig. El (b)). Thus, it was clear that the sloshing modes grew and that no spiral modes 
originated with numerical errors emerged. 

Figure [7] presents the Fourier-transformed normalized amplitudes |c™(c<j)| of sloshing 
m = modes for 2D Models BO, Bl, and B2. The m > modes are plotted on the right side 
of the graph, and the m < modes are plotted on the left side. In this figure, the plots on 
both sides are same due to m = 0, in other words, the left side is just a copy of the right side. 
The results of the new mode analysis indicated that the lower modes tended to be dominant, 
consistent with the feature obtained by the previous mode analysis presented in Fig. H] (a). 
The magnitude of the maximum amplitudes of I = 1 mode and its oscillation frequency uj 
were almost the same between non-rotational Model B0 and rotational Models Bl and B2. 
These results indicated that the rotation hardly affected the characteristics of the sloshing 
modes in SASI. However, the steady I = 2,m = mode was much smaller for rotational 
Models Bl and B2 than for non-rotational Model B0. This may come from the fact that 
centrifugal force acted more strongly on the matter around the equatorial plane than near the 
poles. The vertically long sllipsoidal shock wave oscillating above and below for Model B0 
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might be transformed into the more spherical configuration of Models Bl and B2 by rotation. 



3.2. Non-Axisymmetric Rotational Model 

Next, we considered the results of non-axisymmetric flows. Here, we call the non- 
axisymmetric flow as the 3D flow. When we imposed a random perturbation on the non- 
rotational flow, both axisymmetric m = modes and non-axisymmetric m/fl modes grew 



(llwakami et al.l 120081 ) . In this situation, both sloshing modes and spiral modes could be 
generated. However, our previous mode analysis was not able to distinguish between slosh- 
ing modes and spiral modes for non-axisymmetric m^O modes, because the spectra were 
obtained from the instantaneous deformation of the shock wave. So we used a new mode 
analysis in an attempt to confirm whether spiral modes were actually generated. Further- 
more, we investigated the effect of rotation on sloshing and spiral modes in the nonlinear 
phase. 

Figure [8] plots the average shock radius Rs as a function of time for 3D Models CO, CI, 
and C2. As with 2D Model BO, the average shock radius increased with growing SASI from 
100ms, and then it remained nearly constant until 500ms. At 500ms, the flow started to 
come under the influence of rotational flow falling to PNS for Models CI and C2. Unlike the 
2D models, the faster the rotation, the larger the average shock radius tended to be. When 
we imposed more rapid rotation to the flow than in Model C2, the expansion of the shock 
wave continuously proceeded to the outer boundary. In this study, we therefore focused on 
only the non-exploding models of Models CI and C2. 

Figure M illustrates the flow fields for non-rotational 3D Model CO and for rotational 
3D Model CI. Symmetry with respect to the z-axis was broken for both models (Figs. [9](a), 
(c)). Many large high-entropy blobs were observed inside the shock wave. The blobs with 
circulating flow arised and expanded with growing SASI, and pushed away the shock wave 
outward. Matter infalling from a triple point on the shock wave ran through the interstices of 
the blobs with high velocity, and the stream accreted on the PNS in an arc (Figs. E](b), (d)). 
A triple point is the connection point (or segment) of two shock waves propagating from both 
sides. When rotation was not imposed, no specific rotation axis was observed in the flow. On 
the other hand, as a result of introducing the rotation whose axis corresponded to the z-axis, 
the flow inside the shock wave rotated globally around the z-axis, and accreting spiral flows 
had higher velocity than that for non-rotational models because of <p component of velocity. 
After the accreting flows turned in the rotational direction near the PNS, the flows went 
up inside the high-entropy blobs with higher velocity than that for non-rotational models 
(Figs. M (d))- This is favorable for the shock expansion. Moreover, after the upward flow with 
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high velocity turned in the rotational direction along the shock front, the centrifugal forces 
acted on the flow so that the shock wave could be pushed away furthermore (Figs. 0(d)). 
These flow behaviors might be why the average shock radius for rotational 3D models tended 
to be larger than that for non-rotational ones (Fig. [S]), while the average shock radii were 
not different so much between rotational and non-rotational 2D models (Fig. [5]). The flows 
accreting from the triple points with high velocity were developed in the 9 direction for 2D 
models (Fig. [6] (a)). Thus, the flows running in the rotational direction for 2D rotational 
models do not have so much high velocity as that for 3D ones (Fig. 0(b)). Centrifugal forces, 
therefore, less acted on the 2D flow than the 3D one so that the shock wave might hardly 
expanded for 2D models. 

Figure [TU1 depicts the Fourier-transformed normalized amplitudes |c™(a;)| of sloshing 
m = modes. Figure [TD] (a) presents the results of non-rotational 3D Model CO. Lower 
modes were dominant. The amplitudes for 3D Model CO tended to be smaller than those 
for 2D Model BO (Figs. [7] (a) and Fig. [10] (a)). These results also did not conflict with those 
of previous mode analyses (Figs. 0] (a), (b)). Figures [TUl (b) and (c) present the results of 
rotational 3D Models CI and C2. The steady I — 2, m — sloshing modes grew significantly. 
The spherical shock wave for non-rotational models might change into a horizontally long 
ellipsoidal one because of the centrifugal forces. The amplitudes of the other sloshing modes 
also became somewhat large in response to the increase of rotational velocity, perhaps due 
to the nonlinear effects of coupling the m = and m/0 modes. 

Figure [TT] demonstrates the Fourier-transformed normalized amplitudes |c™(u;)| of \m\ — 
I modes for 3D Models CO, CI, and C2. The basis function of \m\ = I modes did not have 
any node point in the 9 direction. The amplitudes of the m > modes are plotted on 
the right side of the graph, and those of the m < modes are plotted on the left side. If 
the distribution of |c™(u;)| is completely symmetric with respect to uo = in this figure, 
this result can be interpreted in three ways: (1) only sloshing ±m modes (i.e., degenerate 
spiral ±m modes) grow; (2) a spiral +m mode and a spiral — m mode emerge with the 
same probability; (3) both (1) and (2) occur during sampling time. On the other hand, an 
asymmetric distribution of |cj"(w)| suggests that the spiral m modes with a larger amplitude 
is dominant. Figure [TT] (a) plots the non-rotational 3D Model CO, revealing asymmetries in 
the distribution. We recognized that spiral modes were generated in the randomly perturbed 
flow even for the non-rotational model. However, the maximum amplitudes were roughly 
same between m > modes and m < modes. We assumed that sloshing ±m modes 
dominated, or spiral +m modes and spiral — m modes would appear almost equally from 
600ms to 1000 ms. Figures [TT] (b) and (c) present the results of rotational Models CI and 
C2. The m < modes grew more than the m > modes. Furthermore, the amplitude of the 
m < modes became larger with increasing rotation. We introduced the clockwise rotation 
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in the same direction as m < modes. Therefore, the faster the rotation, the larger the 
amplitudes of the mode rotating in the same direction as a globally rotational flow tends to 
be. 



4. Summary and Discussion 

We investigated the effects of rotation on SASI in the nonlinear phase with three- 
dimensional hydrodynamics simulations, for the purpose of application to the supernova 
core in the postbounce phase. When the SASI entered the nonlinear phase, we imposed 
rigid rotation at the outer boundary of the iron core, whose specific angular momentum 
on the equatorial plane was assumed to agree with the recent stellar evolution calculations 
with magnetic fields. After the rotational flow arrived at the shock wave, rotation began 
to influence the SASI in the nonlinear phase. Focusing on this stage, we performed mode 
analysis for the nonspherical deformation of the shock front, using spherical harmonics in 
space and Fourier decompositions in time. 

First, we examined the effects of rotation on SASI for the axisymmetric flow in which 
only sloshing modes existed before rotation was added. In the ranges of rotational strength 
enough not to explode, the average shock radius hardly changed with increasing angular 
momentum. Moreover, mode analysis revealed that the sloshing modes were also insensitive 
to rotation except for the steady / = 2, m = mode. Co mbining our results w i th th e 



outcomes obtained in the linear analyses by lLamingi (120071 ) and lYamasaki fc Foglizzol (120081 ) 



rotation should barely affect the growth rates of sloshing modes. 

Moreover, we studied the effects of rotation on SASI for the non-axisymmetric flow in 
which both sloshing and spiral modes existed before rotation was added. In contrast with 
the axisymmetric models, the shock radius tended to expand more with increased rotation. 
As faster rotation was added to the flows, spiral flows ran with higher velocity from the triple 
points to near the PNS, and then flows circulated with higher velocity inside the high-entropy 
blobs. Large centrifugal force acted on the circulating flows inside the blobs, and the flows 
running in the rotational direction might push the shock wave further outward. As with 
the axisymmetric models, the effect of rotation on the sloshing modes became prominent 
for the steady I = 2, m = mode. We observed that the other sloshing modes also grew 
slightly with increasing rotational velocity. This result might be due to nonlinear effects of 
the coupling of the m ^ modes and m = modes. For non-rotational models, sloshing 
±m modes were dominant, or spiral +m modes and spiral — m modes emerged with almost 
the same probability. However, for rotational models, the spiral modes rotating in the same 
direction as the rotational flow developed significantly with faster rotation. These results 
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agreed with the linear analyses of SASI by [Laming (120071 ) and lYamasaki Foglizzol (120081 ). 



It should be noted that the simulations highlighted in this paper are a first step towards 
realistic three-dimensional modeling of supernova explosions. The approximations adopted 
in this paper (e.g., the replacement of the PNS by the fixed inner boundary and the light- 
bulb approach with the constant neutrino luminosity) need improvement. It is important 
to clarify how rotational SASI impacts neutrino heating, but transport schemes beyond the 
light-bulb approaches are unquestionably needed. Confirming the outcomes of this paper 
will require consistent simulations in 3D, covering the entire stellar core and starting from 
gravitational collapse with better neutrino transport, which is computationally prohibitive at 
present. Additionally, understanding the effects of rotation on the linear growth of SASI and 
constructing steady accretion flows with rota tion are major undertaki n gs. T he generation of 



pulsar spins by SASI has been addressed by iBlondin fc Mezzacappal (120071 ). Currently, we 



are systematically investigating a possible correlation between the kick and spin of PNS, and 
our results will be presented in a forthcoming paper, along with the discussions of magnetic 
effect on SASI (Iwakami et al. in preparation). 
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Table 1: Summary of all models. 



Model Perturbation (3% L [cm 2 s 1 ] b 



AO 


none 








Al 


none 




0.0100 


4.6 x 10 15 


A2 


none 




0.0125 


5.8 x 10 15 


A3 


none 




0.0150 


6.9 x 10 15 


A4 


none 




0.0200 


9.2 x 10 15 


BO 


1=1,171= 









Bl 


1=1,171= 





0.0100 


4.6 x 10 15 


B2 


1=1,171= 





0.0125 


5.8 x 10 15 


CO 


random 








CI 


random 




0.0100 


4.6 x 10 15 


C2 


random 




0.0150 


6.9 x 10 15 



"Parameter for rotation in Eq. 2] 

& Specific angular momentum on the equatorial plane. 
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0.2 




Fig. 1. — Radial distributions of the component of velocity along the equatorial plane 
(6 = 90°, = 0°) for Model Al. The initial flow was the steady spherically symmetric one. 
The rotation for (3$ = 0.01 is imposed on the outer boundary at t — 0ms. The plots of 
t = 200ms and t = 300ms overlap. 
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Fig. 2. — Time evolutions of the average shock radius Rs in the unperturbed flow for Models 
AO, Al, A2, A3, and A4. The time denoted by a dashed line corresponds to the time of 
rotation arrival at the inner boundary. 
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(a) Model AO {b}MtnidAi 



Fig. 3. — Iso-entropy surfaces and velocity vectors in the equatorial section at t = 300ms in 
the unperturbed flow (a) without rotation for Model AO and (b) with rotation for Model Al. 
The length of each side of a panel corresponds to 5.0 x 10 7 cm. 
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t [ms] 

(a) Model Bl 




Z 200 400 600 800 1000 

t [ms] 

(b) Model CI 



Fig. 4.— Time evolutions of the normalized amplitudes |c™(t)/c[J()f:)|. We imposed (a) 
axisymmetric I — 1, m = perturbation for 2D Model Bl and (b) non-axisymmetric random 
perturbation for 3D Model CI on the spherical symmetric flow at t — 0ms. The dotted 
line indicates the starting time to impose rotation on the bouter boundary, the dashed line 
indicates the time of rotation arrival at the inner boundary, and the dot-dashed line indicates 
the onset of the rotation phase. 




Fig. 5. — Time evolutions of the average shock radius Rs without rotation for Model BO 
and with rotation for Models Bl and B2. The times denoted by the dotted line, dashed line, 
and dotted-dashed line are the same as in Fig. HI 
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Fig. 6. — Partial cutaway view of the iso-entropy surfaces and the velocity vectors on the 
cutting plane at t — 917ms with rotation for Model Bl. (a) The object having three cutting 
planes (9 = 90°, = 0°,250°) viewed from the side of — y direction and (b) its equatorial 
section of 9 = 90° viewed from z direction. 
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Fig. 7.— Fourier-transformed normalized amplitude |c™(c<j)| of sloshing m = modes 
without rotation for (a) Model B0 and with rotation for (b) Model Bl and (c) Model B2. 
Normalized amplitudes were estimated by integrating the shock deformation during the 
rotation phase. 
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Fig. 8. — Time evolutions of the average shock radius Rs without rotation for Model CO 
and with rotation for Models CI and C2. 



-25 - 




iclMvdcK'l fd.1 Model CI 



Fig. 9.— Partial cutaway view of the iso-entropy surfaces and the velocity vectors on 
the cutting plane at t — 800ms (a), (b) without rotation for Models CO and (c), (d) with 
rotation for Model CI. One can see (a), (c) the object having three cutting planes from the 
— y direction and (b), (d) its equatorial section from z direction. 
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Fig. 10. — Fourier-transformed normalized amplitude |c™(co>)| of sloshing m = modes 
without rotation for (a) Model CO and with rotation for (b) Model CI and (c) Model C2. 
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Fig. 11. — Fourier-transformed normalized amplitude |c["(u;)| of \m\ = I modes without 
rotation for (a) Model CO and with rotation for (b) Model CI and (c) Model C2. 



